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Abstract 

We study a genealogical model for continuous-state branching processes with immigration 
with a (sub)critical branching mechanism. This model allows the immigrants to be on the 
same line of descent. The corresponding family tree is an ordered rooted continuum random 
tree with a single infinite end defined by two continuous processes denoted by (H t ;t > 0) 
and (H t ;t > 0) that code the parts at resp. the left and the right hand side of the infinite 
line of descent of the tree. These processes are called the left and the right height processes. 
We define their local time processes via an approximation procedure and we prove that they 
enjoy a Ray-Knight property. We also discuss the important special case corresponding to the 
size-biased Galton- Watson tree in the continuous setting. In the last part of the paper we give 
a convergence result under general assumptions for rescaled discrete left and right contour 
processes of sequences of Galton- Watson trees with immigration. We also provide a strong 
invariance principle for a sequence of rescaled Galton- Watson processes with immigration 
that also holds in the supercritical case. 
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1 Introduction 

1.1 The genealogy of Galton- Watson branching processes. 

The continuous-state branching processes with immigration (CSBPI for short) have been intro- 
duced by Kawasu and Watanabe in [19]. They are the continuous analogues of the Galton- Watson 
processes with immigration. In this paper we discuss a genealogical model for CSBPI's that can 
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be described in the discrete setting as follows: let /i and v be two probability measures on the 
set of non-negative integers denoted by N. Recall that a Galton- Watson process Z = (Z n ; n > 0) 
with offspring distribution fi and immigration distribution v (a GWI(/x, i/)-process for short) is 
an N-valued Markov chain whose transition probabilities are characterized by 

E [x z ^ | Z n ] = g(x) z "f (x) , x G [0, 1], (1) 

where g (resp. /) stands for the generating function of (i (resp. v). 

The genealogical model we consider can be informally described as follows. Consider a pop- 
ulation evolving at random roughly speaking according to a GWI(/J, ^)-process: The population 
can be decomposed in two kinds of individuals, namely the mutants and the non-mutants; there 
is exactly one mutant at a given generation; each individual gives birth to an independent num- 
ber of children: the mutants in accordance with v and the non-mutants in accordance with \x. 
We require that all the mutants are on the same infinite line of descent. Except in one part of 
the paper we restrict our attention to a critical or subcritical offspring distribution fi: 

£ < 1. 

fc>0 

Then the resulting family tree is a tree with a single infinite end. We call such trees sin-trees 
following Aldous' terminology in [2] in Section 4. In order to code them by real-valued functions, 
all the discrete trees that we consider are ordered and rooted or equivalently are planar graphs 
(see Section 2.1 for more details). So we have to specify for each mutant how many of its children 
are on the left hand of the infinite line of descent and how many are on the right hand. We 
choose to dispatch them independently at random in accordance with a probability measure r 
on the set {(k,j) € N* x N* : 1 < j < k}. More precisely, with probability r(k,j) a mutant has 
k — 1 non-mutant children and one mutant child who is the j-th child to be born; consequently 
there is j — 1 non-mutants children on the left hand of the infinite line of descent and k — j on 
the right hand. The immigration distribution v then is given by 

u(k-i) = r ( fc >-?) . k ^ L 

l<j<k 

The probability measure r is called the dispatching distribution and the resulting random tree is 
called a (p,r)- Galton- Wats on tree with immigration (a GWI(/x, r)-tree for short). Indeed, if we 
denote by Z n the number of non-mutants at generation n in the tree, then it is easy to see that 
Z = (Z n ; n > 0) is GWI(/U, i/)-process. 

Let us mention the special case r(k,j) = fi(k)/fl, 1 < j < k that corresponds to the size- 
biased Galton- Watson tree with offspring distribution /j,. This random tree arises naturally by 
conditioning (sub)critical GW-trees on non-extinction: see [1], [2], [20], [21], [28] and in the 
continuous case [26]. 

We shall code a GWI(/x, r)-tree r by two real-valued functions in the following way: think 
of r as a planar graph embedded in the clockwise oriented half-plane with unit edge length and 
consider a particle visiting continuously the edges of r at speed one from the left to the right, 
going backward as less as possible; we denote by C s (r) the distance from the root of the particle 
at time s and we call the resulting process C (r) := (C s (r); s > 0) the left contour process of r. 
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It is clear that the particle never reaches the right part of r but C (r) completely codes the left 
part of r. We denote by C (r) the process corresponding to a particle visiting r from the right to 
the left so we can then reconstruct r from (C (r), C (r)): see Section 2.2 for precise definitions 
and other codings of sin-trees. 



1.2 Background on continuous-states branching processes with immigration. 



The main purpose of the paper is to provide a genealogical model for CSBPI's and to build a con- 
tinuous family tree coded by two functions playing the role of C and C . Before discussing it more 
specifically, let us recall from [19] that continuous-state branching processes with immigration 
are [0, oo]-valued stochastically continuous Markov processes whose distribution is characterized 
by two functions on [0, oo): a branching mechanism ip such that (— ip) is the Laplace exponent of 
a spectrally positive Levy process denoted by X = (X t ;t > 0) and an immigration mechanism 
ip that is the Laplace exponent of a subordinator denoted by W = (Wt]t > 0): 



E 



-Alt 



E 



-xw t 



= e -MA) 



X,t > 0. 



More precisely, Y* = (Y t *;t > 0) is a (ip, (^-continuous-state branching process with immigration 
(a CSBPI(ip, ip) for short) if its transition kernels are characterized by 



E 



e - XY a + b \Y* 



exp I —u(a, \)Yfr 



ds <p(u(s, A)) I , a, b, A > 0, 



where u(a, A) is the unique nonnegative solution of the differential equation 
d 

—u(a, A) = — ip (u(a, A)) and u(0, A) = A , a, A > 0. 



Note that this differential equation is equivalent to the integral equation 

rX du 

lu(a,\) 



/ 

Jut 



a. 



(2) 



(3) 



(4) 



Observe that oo is an absorbing state. In the paper, we will only consider conservative processes, 
that is processes such that a.s. Y t * < oo, t > 0. This is equivalent to the analytical conditions 



du 



o+ 



\*P(u)\ 



oo and ^(0) = ip(0) = . 



(5) 



Observe that if there is no immigration, that is if <p = 0, the process is simply a V'-continuous- 
state branching process (a CSBP(^) for short). We shall denote CSBPI's and CSBP's in a 
generic way by resp. Y* = (Y t *;t > 0) and Y = (Y t ;t > 0). We refer to [19] or [29] for a precise 
discussion of CSBPI's and to [7] and [15] for results on CSBP's. 

Except in Section 4.1 we only consider CSBPI's with (sub)critical branching mechanism. 
This assumption is equivalent to the fact that X does not drift to +oo. In that case ip is of the 
form 

iP(\) = a\ + l3\ 2 + f ir(dr)(e~ Xr - 1 + Ar) , A > 0, (6) 

J(0,oo) 
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where a, (3 > and 7r is a cr-finite measure on (0, oo) such that Jj ^ ir(dr)(r A r 2 ) < oo. We 
also assume 



which is equivalent to the a.s. extinction of CSBP(?/>) (see [15] for details). Note that (7) implies 
at least one of the two conditions: [3 > or Jq rir(dr) = oo that guarantee that X has infinite 
variation sample paths (see [5] for details). 

We build the family tree corresponding to a CSBPI(^, (p) thanks to two continuous processes 
(H t ;t > 0) and (H t ;t > 0) called the left and the right height processes that are viewed as the 
contour processes of the parts at the left and at the right hand of the infinite line of descent. 
More precisely, our construction relies on two auxilliary processes: 

• The first one is the height process H = (H t ;t > 0) introduced by Le Gall and Le Jan [25] 
coding the genealogy of ^-continuous-state branching processes (see also [11] for related results). 
H is obtained as a functional H(X) of the spectrally positive Levy process X with exponent 
V>. More precisely, for every t > 0, H t "measures" in a local time sense the size of the set 
{s < t : X s _ = inf[ s t ] X r } (see Section 3.1 for a more precise definition). Assumption (7) 
is equivalent for H to have a continuous modification. From now on, we only consider this 
modification. 

An important role is played by the excursion measure N of X above its minimum process. In 
the quadratic branching case ^{u) = cm 2 , X is a (scaled) Brownian motion, the height process 
H is a reflected Brownian motion and the "law" of H under N is the Ito measure of positive 
excursions of the linear Brownian motion: This is related to the fact that the contour process of 
Aldous' Continuum Random Tree is given by a normalized Brownian excursion (see [3] and [4]), 
or to the Brownian snake construction of superprocesses with quadratic branching mechanism 
(see e.g. [24]). For a general tp, the limit theorems for the contour processes of discrete Galton- 
Watson trees given in [11], Chapter 2 and the Ray-Knight property for the local times of H 
proved in Theorem 1.4.1 [11] both strongly justify that the height process is the right object to 
code the genealogy of (sub)critical CSBP's and that H under the excursion measure N is the 
contour process of a continuum random tree that is called the Levy tree: we refer to [12] for a 
precise definition of Levy trees in term of random metric M-trees space (see also [13] for related 
topics). All the results about height processes used in the paper are recalled in Section 3.1. 

• The second process is a bivariate subordinator (U, V) = ((Ut,Vt); t > 0), namely a [0, oo) x 
[0, oo)-valued Levy process started at 0. Its distribution is characterized by its Laplace exponent 
$: 



Levy measure R plays the role of the discrete dispatching measure r. Roughly speaking, think of 
the population as being indexed by the positive real numbers and let us make an informal analogy 
with the discrete model: If the height t G [0, oo) in the family tree of the CSBPI corresponds to 




(7) 



E [exp (-pU t - qV t )] = exp(-t$(p, q)). 



<3? is of the form 
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generation n in the discrete GWI-tree, then Ut (resp. Vt) corresponds to the sum of the numbers 
of immigrants at the left (resp. the right) hand of the infinite line of descent from generation 
to generation n. Then, a jump of (U, V) occuring at time t corresponds to a total amount 
Ut — Ut- + Vt — Vt- of immigrants arriving at height t in the family tree: Ut — Ut- of them are put 
at the left hand of the infinite line of descent and they are the initial population of a CSBP(V>) 
; Vt — Vt- of them are put at the right hand of the infinite line of descent and they are also 
the initial population of an independent CSBP (?/>). It implies that the real- valued subordinator 
U + V has Laplace exponent <p and thus, 

<f(p) = $(P,P)- 

More precisely, we define H and H as follows. Let us first introduce the right-continuous 
inverses of U and V: 

Uf 1 = inf{s > : U s > t} and Vf 1 = inf{s > : V s > t}, 

with the convention inf = oo. Let H be the height process associated with the Levy process 
X with Laplace exponent ip and let (H',X') be an independent copy of (H,X). We set for any 
t > 0, It = inf[o,t] X and I[ = inf[ 0) t] X' . Then, we define H and H by 

H t = H t + UZI and H t = H[ + Vz], , t > 0. (8) 

The processes H and H are called respectively the left and the right height processes. The left 
and right height processes are continuous iff U~ l and V^ 1 are continuous, which happens iff U 
and V are not Poisson processes. This is equivalent to the analytical condition 

or i?((0,oo) 2 ) = oo. (9) 

As for Levy trees, it is possible to build a {i^, ^-immigration Levy tree via the left and right 
height processes: To each s£R corresponds a vertex in the continuum tree at height 

J s = l(_ OOi0 )(s) J ff- s + l[ ,oo)(s)H s . 

Suppose that s < s' . The common ancestor of the vertices corresponding to s and s' is situated 
at height m(s,s') = inf{J u ; u € I(s,s')}, where I(s,s') is taken as [s,s'] if ^ [s,s'] and as 
R\ [s, s'] otherwise. Then, the distance separating the vertices corresponding to s and s' is given 
by 

d(s, s') = J s + J s i — 2m(s, s'). 

Check that d is a pseudo-metric on R. We say that two real numbers s and s' are equivalent if 
they correpond to the same vertex, that is: d(s, s') = 0. This equivalence relation is denoted by 
s ~ s' and we formally define the (tp, ^-immigration Levy tree as the quotient set T* = R/ ~ 
equipped with the metric d that makes it a random Polish space. Arguing as in [12], we can 
show that a.s. (T*,d) is a real tree. 

This genealogical model is clearly related to the model discussed in [22] by A. Lambert where 
all the population is on the left hand side of the infinite line of descent so that only one contour 
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process is needed to encode the family tree of the CSBPI's. In Lambert's paper, this contour 
process is defined as a functional of a Markov process X* generalizing the functional giving the 
height process introduced by Le Gall and Le Jan. This Markov process can be constructed either 
pathwise in terms of a spectrally positive Levy process and an independent subordinator, or in 
distribution thanks to Ito's synthesis theorem. A. Lambert also defines in a weak sense the local 
time processes of the resulting contour process and states a generalized Ray-Knight theorem by 
proving they are distributed as a CSBPI. 

1.3 Statements of the main results. 

In this paper, the model that we consider allows to have population on both sides of the infinite 
line of descent which turns out to be a natural case to discuss for we can define continuous 
analogues of the discrete size-biased trees. In particular we show in Theorem 1.3 stated below 
that, as in the discrete case, the continuous analogues of size-biased trees are the limit of the 
family trees of (sub) critical CSBP's conditioned on non-extinction. This special case strongly 
motivated the present work. In the more general case, we also provide strong approximation for 
the local time processes of the right and the left contour processes (Proposition 1.1) and we state 
the Ray-Knight property for them (Theorem 1.2). We prove functional weak convergence for the 
rescaled contour processes of discrete GWI-trees (Theorem 1.5). We also state in Theorem 1.4 a 
general limit theorem for GWI processes (with possibly supercritical offspring distribution), in 
the same vein as Theorem 3.4 [17] of Grimvall. 

More specifically, let us now give a detailed presentation of the main results of the paper. The 
first one defines the local time processes of the two contour processes by a strong approximation 
procedure. 

Proposition 1.1 Assume that (7) and (9) hold. Then, there exist two jointly measurable pro- 
cesses { L a s ; a,s>0) and (L^; a, s > 0) such that: 

(i) A.s. for any a > 0, (L®; s > 0) and (L®; s > 0) are continuous nondecreasing processes. 



(ii) For all T > the following limit holds in probability: 



sup 

t<T 



0. 



The same limit holds for H and L a . 
(Hi) A.s. for any continuous function g on [0,oo) with compact support and for any t > 0, 

rt ^ roo ^ f't ^ /*oo ^ 

/ dsg(H s )= I daLlg(a) and / dsg(H s )= / daL^g(a). 
Jo Jo Jo Jo 

We next show that the local time processes of H and H enjoy a "Ray-Knight" property: As U 
and V both drift to infinity, so do the left and the right height processes and it makes sense to 
define 

= lim lZ a s and = lim ~L" , a > 0. 

s— +oo s— >oo 

Then, the Ray-Knight theorem can be stated as follows. 
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Theorem 1.2 Assume that (7) and (9) hold. Then, the process (L^+LJ^,; a > 0) is a 
CSBPI(^,if) started at with tp(X) = $(A, A) , A > 0. 

Proposition 1.1 and Theorem 1.2 are proved in Section 3.2 while Section 3.3 is devoted to the 
study of the continuous analogue of size biased GW-trees, that corresponds to (tp, ^-immigration 
Levy trees where $ is given by 

p-q 

Here, we have set ^*(A) = ^(A) — a A and when p = q, the ratio (ip*(p) — ~4>*(q))/(p — q) should 
be interpreted as ip'(p) — a. So U + V is a subordinator with Laplace exponent <p = tp' —a that is 
the immigration mechanism of the underlying CSBPI. Now, consider the height process H under 
its excursion measure N and denote by £ the duration of the excursion. As a consequence of (7) 
we get A^(sup sg [o j( j] H s > a) £ (0, oo) for any a > 0. Thus, we can define the probability measure 
iV( a ) = N( ■ | sup H > a) (see Section 3.1 for details). The main result proved in Section 3.3 can 
be stated as follows. 

Theorem 1.3 Assume that (7) holds. Then, 

(H tA( , H {c _ t)+ ; t > 0) under N (a) —— » (H t ,H t ;t> o) 
weakly in C([0, oo), IR 2 ). 

(Here (x) + stands for the non-negative part max(0, x) of x.) The proof of this theorem relies on 
a lemma (Lemma 3.2) that is stated in Section 3.3 and that is an easy consequence of Lemma 
3.4 in [12]. Let us mention that Lemma 3.2 is a generalization of Bismut's decomposition of the 
Brownian excursion. 

Section 4.1 and 4.2 are devoted to limit theorems for rescaled GWI processes and contours 
of GWI trees. The main result proved in Section 4.1 is a strong invariance principle for GWI 
processes. In Section 4.1 and only in Section 4.1 we do not anymore restrict our attention 
to (sub)critical GWI-processes. More precisely, let {ji p ;p > 1) and {y p \p > 1) be any sequences 
of probability measures on N and let x £ [0, oo). We denote by (Yn P ; n > 0), p > 1 a sequence of 
GWI(/i p , f p )-processes started at Yq' p = [px] and we denote by {j p ;p > 1) an increasing sequence 
of positive integers. 

Theorem 1.4 The three following assertions are equivalent: 
(i) For any t > the following convergence 

'-na^- (io) 

holds in distribution in R; Here (Z£ ; t > 0) stands for a non-constant and stochastically 
continuous process such that 

Vt > , P(Z* > 0) > and P(Z* < oo) = 1. 
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(ii) We can find a non-constant spectrally positive Levy process X = (Xt;t > 0) with exponent 
tp and a subordinator W = {Wt',t > 0) with exponent <p such that (5) holds and such that the 
following convergences 

/%> ^ P(Xi G • ) and v p ( -V" ^ P(W 1 € • ) (11) 

hold in distribution in M. (here * denotes the convolution product of measures). 

(Hi) There exists a non-constant and conservative CSBPI(ip,ip) denoted by Y* = (Y t *;t > 0), 

started at Yq = x and such that 

weakly in the cadlag functions space D([0, oo),M) endowed with the Skorohod topology. 

In regard of Theorem 3.4 [17] due to Grimvall that concerns limits of GW-processes without 
immigration, the latter limit theorem is very natural. However it turns out to be new. To 
prove (ii) (Hi) we adapt an argument contained in the proof of Theorem 3.4 [17]; our main 
contribution is the proof of (i) => (ii). 

In Section 4.2 we prove a limit theorem for the genealogy of a sequence of (sub)critical GWI- 
processes: let (n p ;p > 1) be a sequence of offspring distributions such that fi p = Xlfc>o ^Mp(^) < 1 
and denote by g^ the corresponding generating functions. Define recursively gffl by g^ = 
9n-i ° 9^ with g^ = Id. Let (r p ;p > 1) be a sequence of dispatching distributions and let r p 
be a G WI(/x p , r p )-tree. For any n > 0, we also denote by Y£ ' p the number of non-mutants at 
generation n in t p . Recall that (7 P ;p > 1) stands for an increasing sequence of positive integers. 
We suppose that 

li P (—) ^ ^ P(Xi € ■ ) and r p (- , -V" € • ; V 1 € • ) (13) 

\ P J P^°° \p p J p^oo 

weakly in resp. R and M 2 . Here X stands for a spectrally positive Levy process whose exponent 
ip satisfies (6) and (7), and (U, V) stands for a bivariate subordinator whose exponent $ satisfies 
(9). We also make the additional assumption 

which implies that extinction times of GW(// p )-processes converge in distribution in the 7" 1 
time-scale. (14) turns out to be a necessary condition in order to have a strong convergence 
of contour processes of rescaled Galton- Watson trees: see [11] Chapter 2 p. 54 for a precise 
discussion of this point. 

Denote by H and H the left and the right height processes associated with a (tp, <3?)- 
immigration Levy tree as defined by (8) and denote by L and L their corresponding local 
times. The main result proved in Section 4.2 is the following. 
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Theorem 1.5 Assume that (13) and (14) hold. For any p > 0, let t p be a GW (n p )-tree. Then, 



in distribution in B([0, 00), R 2 ) xD([0,oo),R). 

This result relies on combinatorial formulas stated in Section 2.2, on Theorem 1.4 and also on 
Theorem 2.3.1 [11] that guarantees a similar convergence for the rescaled contour processes of 
sequences of Galton- Watson trees without immigration. 

The paper is organized as follows: In Section 2.1, we set definitions and notations concerning 
discrete trees. In Section 2.2, we discuss various codings of sin-trees that are used in the proof of 
Theorem 1.5. In Section 3.1 we recall important properties of the height process that are needed 
to prove Proposition 1.1, Theorem 1.2 and Theorem 1.3 in Section 3.2; Section 3.3. Section 4.1 
and Section 4.2 are devoted to the proofs of Theorem 1.4 and Theorem 1.5. 

2 Sin-trees and sin-forests 
2.1 Definitions and examples 

Let N = {0, 1, 2, . . .} be the set of the nonnegative integers, set N* = N \ {0} and denote the set 
of finite words written with positive integers by U := {0} U Un>i(^*) n ■ Let G U be the word 
u\ . . . u n , Ui G N*. We denote the length of u by \u\ : \u\ = n. Let v = v 1 . . . v m G U. Then 
the word uv stands for the concatenation of u and v: uv = u\ . . . u n v± . . . v m . Observe that U is 
totally ordered by the lexicographical order denoted by < . A rooted ordered tree t is a subset 
of U satisfying the following conditions 

(i) G t and is called the root of t. 

(ii) If v G t and if v = uj for j G N*, then, u G t. 

(iii) For every u G t, there exists k u (t) > such that uj G t for every 1 < j < k u (t). 

We denote by T the set of ordered rooted trees. We define on U the genealogical order =<! by 



If u =4 v, we say that u is an ancestor of v. If u is distinct from the root, it has an unique 
predecessor with respect to who is called its parent and who is denoted by u . We define the 
youngest common ancestor of u and v by the ^-maximal element w G U such that w =4 u and 
w =4 v and we denote it by u A v. We also define the distance between u and v by d(u,v) = 
\u\ + \v\ — 2\u A v\ and we use notation [[u,v]\ for the shortest path between u and v. Let t G T 
and u G t. We define the tree t shifted at u by 9 u (t) = {d£D : uv G i] and we denote by [t] u 



p^oo 




Vu, v G U , u =j v ^=^> 3w G U : v = uw. 
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the tree t cut at the node u : [t] u := {u} U{u€t : v Au ^ u}. Observe that [t] u G T. For any 
ui, . . . , u k G t we also set [t] Uu ... iUk := Mm fl . . . n [i] Ufc and 

[An = [t]{uet: \u\=n} = {u £ t : \u\ < n} , n > 0. 

Let us denote by Q the <7-field on T generated by the sets {t G T : u G i} , u G U. All the 
random objects introduced in this paper are defined on an underlying probability space denoted 
by ($7, J 7 , P). A random tree is then a T-Q measurable mapping r : — ► T. We say that a 
sequence of random trees (r^; k > 0) converges in distribution to a random tree r iff 

Vn > 0, Vt G T, P([r k } n = t) ► P([ r ] n = t) 

fc^oo 

and we denote it by r k dlStr > r. 

fc— >oo 

Let jU be a probability distribution on N. We call Galton-Watson tree with offspring distri- 
bution \i (a GW(/x)-tree for short) any F — Q measurable random variable r whose distribution 
is characterized by the following conditions: 

(i) p(Mt) = *) = M») » * > o. 

(ii) For every i > 1 such that fj,(i) ^ 0, the shifted trees #i(t), . . . , 6i(r) under 
P(- | k (r) = i) are independent copies of r under P. 

Remark 2.1 Let u±, G U such that A Uj ^ {ni, . . . ,Uk} , 1 < i,j < A;, and let r be 

a GW(/i)-tree. Then, as a consequence of the definition of GW-trees, conditional on the event 
{ui, . . . ,u k G r}, ui (t), . . . ,9 Uk (r) are i.i.d. GW(/u)-trees independent of [r] ui) ... jttfc . 

We often consider a forest (i.e. a sequence of trees) instead of a single tree. More precisely, 
we define the forest / associated with the sequence of trees (ti ; I > 1) by the set 

/ = {(-l,0)}u|J{(Z,«), «€t,} 
l>i 

and we denote by F the set of forests. Vertex (—1,0) is viewed as a Active root situated at 
generation —1. Let u' = (l,u) G / with I > 1; the height of u' is defined by \u'\ := |tt| and its 
ancestor is defined by (I, 0). For convenience, we denote it by 0/ := (I, 0). As already specified, 
all the ancestors 0i,02, ... are the descendants of (—1,0) and are situated at generation 0. 
Most of the notations concerning trees extend to forests: The lexicographical order < is defined 
on / by taking first the individuals of t\, next those of ti ... etc and leaving (—1, 0) unordered. 
The genealogical order ^ on / is defined tree by tree in an obvious way. Let v' G /. The 
youngest common ancestor of v! and v' is then defined as the ^-maximal element of w' such 
that w' =4 u' and w' ^ v' and we keep denoting it by u' A v'. The number of children of u' is 
k u '(f) := k u (ti) and the forest / shifted at u' is defined as the tree u '(f) '■= G u (ti). We also 
define [f] u > as the forest {u'}L){v' G / : v' Au' / u'} and we extend in an obvious way notations 
[f]u[,...,u' an d [f]n- For convenience of notation, we often identify / with the sequence (ti ; I > 1). 
When (ti ; I > 1) = (t±, . . . , t&, 0, 0, . . .), we say that / is a finite forest with k elements and we 
abusively write / = (ii, . . . , 
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We formally define the set of trees with a single infinite line of descent (called sin-trees for 



T sin = {t G T : Vn > 0, # {v G t : \v\ = n and #6 v (t) = 00} = 1 } . 

Let t G T sin . For any n > 0, we denote by n*(t) the unique individual u on the infinite line 
of descent (i.e. such that #0 u (t) = 00) situated at height n. Observe that u^t) = 0. We 
use notation £oo(t) = {^n(0; ^ > 0} for the infinite line of descent of t and we denote by 
(l n (t) ; n > 1) the sequence of positive integers such that (t) is the word h(t) . . . l n (t) G U. We 
also introduce the set of sin-forests F s j n that is defined as the set of forests / = (ti ; I > 1) such 
that all the trees ti are finite except one sin-tree t; . We extend to sin-forests notations it*, and 
l n by setting l n (f) = l n {t lo ), <(/) = (l ,u n (U ) ) and u^f) = ky 

We now precisely define the Galton- Watson trees with immigration introduced in Section 1: 
Recall that a GWI tree is characterized by 

• its offspring distribution /ionN that we suppose critical or subcritical: p, = Xlfc>o^/ i (^) — 1' 

• its dispatching distribution r defined on the set {(k, j) G N* x N* : 1 < j < k} that prescribes 
the distribution of the number of immigrants and their positions with respect to the infinite line 
of descent. 

More precisely, r is a GWI(^, r)-tree if it satisfies the two following conditions: 

(i) The sequence S = ( (/c M » ( t )(t), Z n+ i(r)) ; n > 0) is i.i.d. with distribution r. 

(ii) Conditional on S, the trees {#n*(r)i( T )> n G N, 1 < i < fe u .( T )(r) with i 7^ l n+ i(r)} 
are mutually independent GW(^)-trees. 

We define a GWI^, r)-forest with I > 1 elements by the forest ip = (t,t\, . . . , ti-±) where the 
Tj's are i.i.d. GW(/u)-trees independent of the GWI(/x, r)-tree r. It will be sometimes convenient 
to insert r at random in the sequence (t±, . . . ,r/_i) but unless otherwise specified we choose to 
put the random sin-tree first in a random sin-forest. 

Example 2.1 : The size-biased GW-tree. Recall from Section 1 that a GW(/i) size-biased tree 
is a GWI(/U, r)-tree with r(k,j) = fi(k)/fi, 1 < j < k. The term "size-biased" can be justified 
by the following elementary result. Let ip be a random forest corresponding to a sequence of I 
independent GW(//)-trees and let <p\, be a GWI(/x, r)-forest with I elements where r is taken as 
above and where the position of the unique random sin-tree in p\, is picked uniformly at random 
among the I possible choices. Check that for any nonnegative measurable functional G on F x U: 



short) by 



E £G ([¥>]„,«) =Y,l» n V[G{lp b ] KM ,u n (p b ))] 



(15) 



_u£tp J n>0 



and in particular 



dP([y b ]» G •) _ Z n {p) 
dP(M„ G •) " ifi n 



where Z n (<p) = #{u G <p : \u\ = n}, n > 0. 
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Example 2.2 : A two-types GW-tree. Let p be probability measure on N x N. Consider 
a population process with two types (say type 1 and type 2) whose branching mechanism is 
described as follows : all the individuals in the tree have the same offspring distribution; namely, 
one individual has k children of type 1 and I children of type 2 with probability p(k, I). We order 
the children putting first those with type 1 and next those with type 2. Assume that we start 
with one ancestor with type 1. If we forget the types, the resulting family tree is a GW(/x)-tree 
where p is given by 

k+l=n 

We assume that p is (sub)critical. For any n > 1, denote by A n the event of a line of descent 
from generation n to the ancestor that only contains individuals with type 1. Then we can prove 
easily 

r under P( • | A n ) d ' str > , 
where Too stands for a GWI(/x, r)-tree where r is given by 

1 k 

r(k,l) = — ^ p(j, k — j) with m = ^ kp(k, N) . 
m j=l k>0 



Example 2.3 : Ascending particle on a GW-tree. Let (ir n ; n > 0) be a sequence of probability 
measures on N such that 7r n ({l, . . . n}) = 1. Let r be a critical or subcritical GW(//)-tree. 
Consider a particle climbing r at random in the following way: it starts at the root 0; suppose 
it is at vertex u € r at time n, then there are two cases: if k u (r) > 0, then at time n + 1 the 
particle goes to v = uj with probability vr fcu ( T )(j); if k u {r) = 0, then the particle stays at u at 
time n + 1. The particle is thus stopped at a final position denoted by U. We can show that 
[r] n conditional on {\U\ > n} is distributed as [roo] n where Too is a GWl(p, r)-tree with 

r(k,l) = ^£) 7^(0, l<l<k. 

Consequently, 

r under P ( • | \U\ > n) d ' str > . 



2.2 Codings of sin-trees 

Let us first recall how to code a finite tree t € T. Let uq = < u\ < . . . < u#t-i be the 
vertices of t listed in the lexicographical order. We define the height process of t by H n (t) = \u n \, 
< n < #t. H{t) clearly characterizes the tree t. 

We also need to code t in a third way by a path D(t) = (D n (t); < n < #i) that is defined 
by D n+ \(t) = D n (t) + k Un (t) — 1 and Do(t) = 0. D{t) is sometimes called the Lukaciewicz path 
associated with t. It is clear that we can reconstruct t from D(t). Observe that the jumps of 
D(t) are not smaller than —1. Moreover D n (t) > for any < n < #t and D# t (t) = — 1. We 
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recall from [25] without proof the following formula that allows to write the height process as a 
functional of D (t): 

H n {t) = # (o < j < n : Dj{t) = inf D k {t)\ , < n < #t. (16) 

L j<k<n J 

Remark 2.2 If r is a critical or subcritical GW(/x)-tree, then it is clear from our definition that 
D(t) is a random walk started at that is stopped at —1 and whose jump distribution is given 
by p(k) = fi(k + 1) , k > —1. Thus (16) allows to write H(t) as a functional of a random walk. 

The previous definition of D and of the height process can be easily extended to a forest 
/ = (ti ; I > 1) of finite trees as follows: Since all the trees t\ are finite, it is possible to list all 
the vertices of / but (—1, 0) in the lexicographical order: uq = 0\ < u\ < . . . by visiting first 
ti, then t2 ... etc. We then simply define the height process of / by H n (f) = \u n \ and D(f) 
by D n+1 (f) = D n {f) + k Un (f) ~ 1 with £>„(/) = 0. Set n p = #t x + . . . + #t p and n = and 
observe that 



H np +k(f) = H k(t P +i) and D np+k (f) = D k (t p+1 ) - p , < k < #t p+1 ,p>0. 

We thus see that the height process of / is the concatenation of the height processes of the trees 
composing /. Moreover the n-th visited vertex u n is in t p iff p = 1 — info<fc<„ D k (f). Then, it is 
easy to check that (16) remains true for every n > when H(t) and D(t) are replaced by resp. 
//(/) and /;(/). 

Let us now explain how to code sin-trees. Let t G T s j„. A particle visiting t in the lexico- 
graphical order never reaches the part of t at the right hand of the infinite line of descent. So 
we need two height processes or equivalently two contour processes to code t. More precisely, 
the left part of t is the set {u G t : 3v € £oo(t) s.t. u < v}. This set can be listed in a lexico- 
graphically increasing sequence of vertices denoted by = uq < u\ < . . . etc. We simply define 
the left height process of t by H n (t) = \u n \ , n > 0. H(t) completely codes the left part of t. 
To code the right part we consider the "mirror image" t* of t. More precisely, let v G t be the 
word c\C2 ■ ■ ■ c n . For any j < n, denote by Vj := c\ . . . Cj the j-th ancestor of v with vo = 0. Set 
c' = k Vj _ x (t) — Cj + 1 and v* = c* . . . c* . We then define t* as {v' , v E t} and we define the right 

height process of t as H(t) := H(t'). 

We next give another way to code a sin tree by two processes called the left contour and the 
right contour processes of the sin-tree t, that are denoted by resp. C (t) and C(t). Informally 
speaking, C (t) is the distance-from-the-root process of a particle starting at the root and moving 
clockwise on t viewed as a unit edge length graph embedded in the oriented half plane. We define 
C (t) as the contour process corresponding to the anti-clockwise journey. So we can also write 
C (t) = C(t'). More precisely, C (t) (resp. C (t)) can be recovered from Hit) (resp. H(t)) 
through the following transform: Set b n = 2n — H n (t) for n > 0. Then observe that 

c ^ = j H n (t) - s + b n if s e [b n , b n+ i - 1), ^ 

1 s - b n+ i + H n+ i{t) if s € [b n+1 - 1, b n+1 ]. 
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The contour process is close to the height process in the following sense: Define a mapping 
q : M + — ► Z + by setting q(s) = n iff s € [b n , b n+ \). Check for every integer m > 1 that 

sup \C s (t)-H q{s) (t)\< sup \C s (t)-H q{s) {t)\ <1+ sup | H n+1 (t)-H n (t)\. (18) 

se[0,m] se[0,6 m ] n<m 

Similarly, it follows from the definition of b n that 

s s i ^ — 

sup |g(a) - -| < sup \q(s) - -| < - sup H n (i) + 1. (19) 

se[0,m] Z se[0,b m ] Z Z n<m 

We now give a decomposition of H(t) and H(t) along £oo(t) that is well suited to G Wl-trees 
and that is used in Section 3.2: Recall that (u n ;n > 0) stands for the sequence of vertices 
of the left part of t listed in the lexicographical order. Let us consider the set {tt*_ 1 (t)i;l < 
i < l n (t);n > 1} of individuals at the left hand of £oo(t) having a brother on £oo(t). To avoid 
trivialities, we assume that this set is not empty and we denote by v\ < vi < ... etc. the 
(possibly finite) sequence of its elements listed in the lexicographical order. 

The forest f(t) = (9 V1 (t),6 V2 (t), . . .) is then composed of the bushes rooted at the left hand of 
£oo(t) listed in the lexicographical order of their roots. Define L n (t) := (li(t) — 1) + . . . + {l n {t) — 1) 
for any n > 1 and Lo(t) = 0; then, consider the p-th individual of f(t) with respect to the 
lexicographical order on f(t); it is easy to check that this individual is in a bush rooted in t at 
height 

a(p) = mf{k > : L k (t) > 1 - inf £>,-(/(*))} . 

j<p 

Thus the corresponding individual in t is u n ( p ) where n(p) is given by 

n(p)=p + a(p) (20) 

(note that the first individual of f(t) is labelled by 0). Conversely, let us consider u n that 
is the n-th individual of the left part of t with respect to the lexicographical order on t. Set 
p( n ) = < n : Uk ^ £oo(t)} that is the number of individuals coming before u n and not 
belonging to £oo(t)- Then 

p(n) = inf{p > : n(p) > n} (21) 
and the desired decomposition follows: 

H n (t) = n-p(n) + H p{n) (f(t)). (22) 

Since n — p(n) = #{0 < k < n : u k € ^oo(^)}, we also get 

a(p(n) - 1) < n - p(n) < a(p(n)). (23) 

Observe that if u n ^ £oo(t), then we actually have n — p(n) = a(p(n)). The proofs of these 
identities follow from simple counting arguments and they are left to the reader . Similar formulas 
hold for H(t) taking P instead of t in (20), (21), (22) and (23). 
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Remark 2.3 The latter decomposition is particularly useful when we consider a GWI(/x, r)- 
tree r: In this case (/(t),/(t*)) is independent of (L(t), L(t')), /(t) and /(r*) are mutually 
independent and /(r) (resp. f(r*)) is a forest of i.i.d. GW(/u)-trees if for a fc > 2 we have 
r(/c, 2) + . . . + r(k, k) ^ (resp. r(fc, fc — 1) + . . . + 1) ^ 0), it is otherwise an empty forest. 
Moreover, the process (L(t), L(t')) is a N x N-valued random walk whose jump distribution is 
given by 

P (L n+ i(r) - L n (r) = m ; L n+1 (r') - L n (r') = m) = r(m + m! + l,m + 1). 

3 Continuum random sin-trees 
3.1 The continuous time height process 

In this section we recall from [25] the definition of the analogue in continuous time of the discrete 
height process defined in Section 2.2. We also recall from [11] several related results used in the 
next sections. 

To define the continuous-time height process, we use an analogue of (16) where the role 
of the random walk is played by a spectrally positive Levy process X = (X t ; t > 0). The 
(sub)criticality of \x corresponds to the fact that X does not drift to +oo. We also assume that 
X has a path of infinite variation (in the finite variation case, the height process is basically a 
discrete process and so is the underlying tree: see [25] and [27] for a discussion with applications 
to queuing processes). As already mentionned in the introduction, this happens if the exponent 
ip of X satisfies conditions (6) and (7). By analogy with (16), the height process H = (H t ; t > 0) 
associated with X is defined in such a way that for every t > H t measures the size of the set: 

{s G [0,t] : AV = inf X r } . (24) 

s<r<t 

To make this precise, we use a time-reversal argument: For any t > 0, we define the Levy process 
reversed at time t by 

xl = X t - X (t _ s) _ , < s < t 

(with the convention Xq_ = 0). Then, X 1 is distributed as X up to time t. Let us set for any 

s > 0, 

S s = sup X r and S\ = sup X\ . 

r<s r<s 

Then , the set (24) is the image of 

{s€[0,t] : St = Xl} 

under the time reversal operation s — > t — s. Recall that under our assumptions S — X is a strong 
Markov process for which is a regular value. So, we can consider its local time process at 
that is denoted by L(X). We define the height process by 

H t = L t {X l ) , t>0. (25) 
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To complete the definition, we still need to specify the normalization of the local time: let us 
introduce the right-continuous inverse of L{X): 



L- 1 = inf{s > : L S (X) > t} 

(with the convention that inf = oo). Define K t by X L -i if t < L^X) and by oo otherwise. A 
classical result of fluctuation theory (see [5] and [6]) asserts that (K t ; t > 0) is a subordinator 
whose Laplace exponent is given by 

E[exp(-XK t )] = exp(-ctV(A)/A) , t,X > 0. 

Here, c is a positive constant that only depends on the normalization of L(X). We fix the 
normalization so that c = 1. When > 0, standard results on subordinators imply for any t > 0, 



H t = im ({5* ; < s < t}) , 



where m stands for the Lebesgue measure on the real line. In particular when X is a Brownian 
motion, we see that H is distributed as a reflected Brownian motion. 

Let us briefly recall the "Ray-Knight theorem" for H (Theorem 1.4.2 [25] and Theorem 1.4.1 
[11]), that can be viewed as a generalization of famous results about linear Brownian motion. 
For any a, t > 0, we introduce the local time of H at time t and at level a that can be defined 
via the following approximation 



lim E 



sup 

Q<s<t 



1 f s 

e Jo 



(26) 



inf{s > : X s = -r} 



r , a > 0. Then, Theorem 1.4.1 [11] asserts that (Y a ; a > 0) is a CSBP(V>) started 



(see Proposition 1.3.3 [11] for details). Next, set for any r > : T, 
and Y a = L% 
at r. 

Although the height process is in general not Markovian, we can still develop an excursion 
theory of H away from 0: Recall notation It = inf s < t X s . Observe that for any t > 0, H t only 
depends on the values taken by X — I on the excursion interval that straddles t. Under our 
assumptions, X — I is a strong Markov process for which is a regular value so that —I can 
be chosen as the local time of X — I at level 0. We denote by iV the corresponding excursion 
measure. Let (gi, d{), i G X be the excursion intervals of X — I above 0. We can check that P-a.s. 

\J{9i, di) = {s>0 : X s - I 8 > 0} = {s > : H s > 0}. 

iei 

Denote by hi(s) = H gi+S , < s < Q = di — gi, i G X the excursions away from 0. Then, each Hi 
can be written as a functional of the corresponding excursion of X — I away from on (gi, di). 
Consequently, the point measure 



M(drdw) = 5^<S(-j fl .,/n)(drdw) 
iei 



(27) 
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is a Poisson point measure with intensity drN(duj). Note that in the Brownian case, N is the 
Ito excursion measure of positive excursions of the reflected linear Brownian motion. 

From now on until the end of this section we argue under N . Let C denote the duration of 
the excursion. The local time processes of the height process {L a s ; 0<s<(),a>0 can be 
defined under iV through the same approximation as before, namely 



lim sup N ly sup 

a>0 V 0<s<tA( 



drl{ a<Hr <a+e} ~ L s 



, t > 0, (28) 



1 

e Jo 

where V is any measurable subset of excursions such that N(V) < oo. The above mentioned 
Ray-Knight Theorem for H implies 

N (l — exp(-AL")) = u(a, A) , a, A > 0, (29) 

where we recall that u is defined by (4). Set v(a) = liniA^oo u(a, A) to be a positive and finite 
quantity by (7) satisfying a = J^ a ^ du/tp(u). By a simple argument discussed in Corollary 
1.4.2[11], we get 

N(L^>0) =N (^supH s > aj = v(a). (30) 

Let a > and set = N( ■ \ sup H > a) that is a well-defined probability measure. The 
Levy tree coded by H under iV7 a ) enjoys a branching property that can be stated as follows: set 



ff = inf{s > : / drl {Hr < a} > t} 
Jo 



We denote by H a the a-field generated by (X^°- , t > 0) and by the class of the iV-negligible sets 
of T . We introduce the excursion intervals of H above a: 

U (gi,d t ) = {s>0 : H s >a} 

ieT(a) 

and we recall from Proposition 1.3.1 [11] the following result. 

Proposition 3.1 (Proposition 1.3.1 [11]) The process (L", s > 0) is measurable with respect to 
TL a . Then, under iV( a ) and conditional on TC a the point measure 

M a {dldu) = 6 (L a g] ,H igt+ . )Ad -a){dldu) (31) 

is independent ofTL a - Moreover it is a Poisson point measure with intensity 

1 [0) l»](0<« N (du). 

Remark 3.1 Proposition 1.3.1 [11] is actually stated under P for the so-called exploration pro- 
cess {pt,t > 0) that is a Markov process taking its values in the space of the finite measures of 
[0, oo) and that is related to the height process in the following way: P-a.s. for any t > 0, the 
topological support of pt is the compact interval [0, H t \. Thus it easy to deduce from Proposition 
1.3.1 [11] a statement for the height process under P and our statement follows from the fact 
that iV( a ) is the distribution under P of the first excursion of H away from that reaches level 
a. 
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3.2 Proof of Theorem 1.2. 

Recall from Section 1 the definition of the left and right height processes H and H of a (if), <&)- 
immigration Levy tree. We first prove Proposition 1.1. 

Proof of Proposition 1.1: We only need to consider H. Recall the notation (gi,di), i £ I 

for the excursion intervals of H away from 0. Set for any a > and any i G I 

(i = di-9i, hi = H gth{ . +di) and a* = (a - UZ] g , )+, 

where for any iflwe have set = x V 0. Recall notation T r = inf{s > : X s = — r} , 

r > 0. Define for any s > 0, 

L J = (-/ s - C/ a _) + AAU a + Y, L a s \ di ~ ■ 

It is clear from the definition that (s,a) —> L™ is jointly measurable. Since the mappings 
s — y {—Is — U a -)+ A AU a and s -> - L^ fl . 

are non-decreasing and continuous, then s — > L " is a non-decreasing mapping and it is continuous 
on every open interval (gi,di) , i G I and also on [Tjj a _, oo). Let s £ [0, Tu a _). Suppose that s 
does not belong to an excursion interval of H away from 0, that is H s = 0. Then, H s = Uz] < a 
and it easy to check that the continuity of H implies the existence of a non-empty open interval 
centered around s on which L a is a constant mapping. These observations imply (i). 

Point (Hi) follows from point (it) by standard arguments. It remains to prove (ii): By (26), 
we see that for any i El, {L^ s+g .^ A ^. — L gi ! a,s>0) only depends on excursion hi. So it makes 

sense to denote it by (L"(/ij) ; a, s > 0). Since P(U a - = U a ) = 1 and by the definition of H we 
a.s. get for any T > 0, 



sup 

0<t<T 



dsl {a<H 3 <a+e} L t 



1 dsl {a . <h .( s) < a .+ e } - Lp(hi 
Jo 



sup 

te[0,TA(i] 



By conditionning on U, we get a.s. 



E 



sup 

0<t<T 



dsl 



{a<H s <a+e} 



T a 



u 



< / dxn e (a- UZ 1 ), (32) 
Jo 



where we have set for any y > 0, 



n t (y) = N l v(y ) sup 

V 0<s<TAC 



e So drl{y<Hr - y+e} ~ LVg ^j 
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with V(y) = {supH > y}. By (30), N(V(y)) = v(y) < oo so (28) applies and we get n e (y) -> 
when e goes to 0, for any fixed y > 0. Moreover, 

n e (y) < N(L y c ) + N (V 1 £ dsl {y<Ha < y+e} ^ = N(L y ( ) + e~ l jT* daN(L a ( ) 



by (28) once again. Then, use (4) and (29) to get 



N(L V _ ^JV (J _ = |-„( M ) |M = e - < 1. 



Thus, n e (y) < 2 and /J 7 " dxn e (a—U x x ) tends a.s. to when e goes to by dominated convergence. 
(ii) follows from (32) by an easy argument. ■ 

Proof of Theorem 1.2: Since a.s. (U, V) has no fixed discontinuity, we can write for any a > 
a.s. 

^c + ^c = E + E L §^) < ( 33 ) 

with an obvious notation for /i', a' and , j € X'. Fix < b\ < . . . < b n . Deduce from (33) 



E 



3 ' J 

"6i i 7 b i 



exp -Ai(L^+ L*)-...-\ n (L%+ L 



OO , 



E 



exp 



u h . 



dxu 



(u- 1 ) - / 

J o 



(34) 
(35) 



where we have set for any s > 0: 

uj(s) = N (l-exp(-\ 1 Lf 1 - sh 



— ... — A n Zi 



(&„—«)+ 



Let (Y" a ; a > 0) denote a CSBP(^) started at Yq = 1. The Ray-Knight property of the local 
times of then implies 



oj(s) = -logE [exp (-AiY" (bl _ s)+ - ... - X n Y {bn _ s)+ )} . 
Now use the Levy-Ito decomposition of (U, V) to get a.s. 



(36) 



f U b„ fVb n 

/ dxcjiU- 1 ) + / dxuiV- 1 ) = 
Jo Jo 

{d + d') j uj{s)ds + E (&U S + AV s )lu(s). 



0<s<b n 



Recall that </?(A) = <E>(A, A) , A > and deduce from the previous identity: 



E 



exp 



(- dxuiU- 1 ) - dxuiV- 1 )^ = exp (- ^{u{s))d. 
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Denote by (Y*; a > 0) a CSBPI(?/> ,<f) started at Yq = 0. An elementary computation (left to 
the reader) shows that 

E [exp (-\ 1 Y* i - ... - X n Y b * n )] = exp (- ^{u{s))ds 

which completes the proof of Theorem 1.2 by (34) and (36). ■ 




Remark 3.2 Let us explain how Theorem 1.2 extends to a CSBPI^,'/ 9 ) started at an arbitrary 
state r > 0. Set H r t = H t + U7} It _ r) . Thus, H° = H. Observe that H r coincides with H up 
to time T r = inf{s > : X s = — r}. Then, use the Markov property at time T r to show that 
^T r +- 1S independent of H r ATr and distributed as H°. Proposition 1.1 and (26) make possible 
to define a local time process for H r denoted by (L r s ' a ; a,s > 0) that satisfies properties (i), 
(ii) and (Hi) of Proposition 1.1. Moreover, the previous observations imply that 

(T^; a > 0) and (Xoo ~ ^ a > 0) 

are two independent processes: the first one is distributed as a CSBP(^) started at r and the 
the second one is a CSBPl(ip,(p) started at 0. Then deduce from (4) and (3) that the sum of 
these two processes is distributed as a CSBPI(^,<^) started at r. 



3.3 Proof of Theorem 1.3. 

In this section we discuss the ^-size-biased Levy tree case, namely 

p-q 

where we have set ip*(X) = V'(A) — a. Let us introduce the last time under level a for the left 
and the right height processes: 

o 7 a = sup{s > : H s < a} and ~a a = sup{s > : H s < a}. 

One important argument in the proof of Theorem 1.3 is the following lemma. 



Lemma 3.2 Assume that (7) holds. Then, for any positive measurable function F and G, 

N (i CdsF ( H -^ G ( H «-)MCs))) = J™ dae- aa E[F(H. A ^ a )G(H.^ a )_ , 
and for any a > 0, 

N^dL:F(H. As )G(H i( _. )A{( _ s) )^=e^ a B [f(H . a ^JG(h . a7?o ) 
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Proof: The second point of the lemma is an easy consequence of the first one and of (28). 
Thus, we only have to prove the first point. To that end, we introduce Mf the space of all finite 
measures on [0, oo). If ji G Mf, we denote by H(fi) € [0, oo] the supremum of the (topological) 
support of fj,. We also introduce a "killing operator" on measures defined as follows. For every 
x > 0, k x fi is the element of Mf such that k x /j,([0, t]) = /x([0, i\) A (/i([0, oo)) — x) + for every 
t > 0. Let Mf stand for the set of all measures \i G Mf such that H{fj) < oo and the topological 
support of n is [0, H(fj,)]. If fx £ Mf, we denote by Q M the law under P of the process defined 
by 

H? = H{k. It p) + H t , if t < T M , 
H? = 0, if t > T M , 

where = inf{i > : X t = — (/x, 1)}. Our assumption \x € Mf guarantees that has 

continuous sample paths, and we can therefore view as a probability measure on the space 
C_|_([0, oo)) of nonnegative continuous functions on [0, oo). For every a > 0, we let M a be the 
probability measure on (Mf) 2 that is the distribution of (l[o, a ](*) dUt, l[o, a ](*) dVt). 

The main argument of the proof is the key-Lemma 3.4 [11] that asserts that for any nonneg- 
ative measurables functions F and G on C+([0, oo)), 



if 



■V( / r/.,F(# (8 _. )+ jG(ff (a+ . )AC 

poo r r 

= J dae- aa J M a (dfidu) J Q^{dh)Q u {dh')F{h)G{h'). (37) 

We use (37) to complete the proof as follows: First observe that (37) implies that the height 
process is reversible under N, namely 

(H 8 ; < s < (l = ] (H c _ s ; < s < Q under N. 
Then fix r > 0. By reversing one-by-one the excursions of H away from on [0, T r ], we get 

(r + /( Tr -.) + - H {Tr _. )+ ) (la ^ ( I-ATr , H. ATr ). (38) 
Next, fix a > and set \i = l[ 0iO ](t) dUt- Note that for any x, 

H(k x n) = a - If- 1 . (39) 
Deduce from (39), (38) and the definition of the left height process that 

(^) rr/U 

A similar identity holds for the right height process and the lemma follows from (37). ■ 

Remark 3.3 This lemma can be viewed as the continuous counterpart of identity (15). 

Remark 3.4 In the Brownian case Y>(A) = A 2 /2, the left and the right height processes are two 
independent three-dimensional Bessel processes and the lemma is a well-known identity due to 
Bismut [8] (See [9] for a generalization to spectrally Levy processes and [10] to general Levy 
processes). 
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Proof of Theorem 1.3: Let b > 0. For any oj in B([0, oo),R) we introduce Tbioj) = inf{s > : 
lo(s) > b} (with the usual convention inf = oo). To simplify notations we set 



H = H, 



(<-•)- 



, n = n(H) and % = T b (H) 



We only have to prove the following convergence for any bounded measurable function F, 



lim N. 



a) (n 



H-/\T b i H-r\T b , 



E 



•At 6 / 



b > 0, 



(40) 



(with an evident notation for and ~Tb) since it implies for any t > 

lim lim N/ a ) (r b , %<t)= lim P (V 6 , ~r b < t) = . 

Let us prove (40): deduce from (28) that N-a.e. the topological support of dL b is included in 
C [t&, C — Tfe]. Thus , by Lemma 3.2 



AT (L\F{H. hTb ,H. h% )) = e~ ab V [f(H 



(41) 



Recall from Section 3.1 the notation (gj,dj), j G X(6) for the excursion intervals of H above level 
b. For any a > b, we set 

Z% = #{jel(b): sup tf s >a} 

s G(9j. rf j) 

that is the number of excursions above level b reaching level a — b. Deduce from Proposition 3.1 
that conditional on Jib under Nr b \ , the random variable Z£ is independent of Tib and distributed 
as a Poisson random variable with parameter L^N(sup H > a — b) = L^v(a — b). Then use (41) 
and the obvious inclusion {sup H > a} C {sup H > b} to get 



N {a) (zSF{H. An ,H. A% j) 



v(a-b) c _ ab 
v(a) 



E 



F(H 



Ati,i H -Kr b ) 



Let C be a bounding constant for F. Then, 



N 



a) 



ZZF(H. ATb ,H. Afb )) - N {a] (F(H. ATb ,H. Afb )) | < CN (a) (Z^, Z a h > 2 



(42) 



(43) 



Now, observe that 



N( a ) (ZS; Z% > 2) = V -^N (L\{\ - exp(-L>v(a - b)))) 
= Vj ^(m<b,0)-lu(b,v(a-b))). 
Since lim^oo v(a — b) = and v(a) = u(b, v(a — 6)), we get 



lim v(a — b)/v(a) = ( — u(b, 0) 



and 



9 



lim — u(b, via — b)) 



d_ 



u(b, 0). 



Thus, linia^oo N^ (Z^; Z£ > 2) = and (40) follows from the latter limits combined with (41), 
(42) and (43). ■ 
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4 Limit theorems 

4.1 Proof of Theorem 1.4. 

Recall the notations of Section 1: Let (fJ> p ;p > 1) and (v p \p > 1) be any sequences of probability 
measures on N. In particular, we do not anymore assume that the /j, p s are (sub)critical. Let 
("f P ',P > 1) be an increasing sequence of positive integers. Denote by and the generating 
functions of resp. fi p and v v . Let x G [0, oo). Recall that for any p > 1, (Yn ,p ; n > 0) stands for 
a GWI(/i p , f p )-process started at Y ' p = [px]. We also need to introduce for any p > 1 a random 
walk (W p ; n > 0) independent of the Y*' p, s, started at and whose jumps distribution is v p . We 
denote by (Y p ;n > 0), p > 1 a sequence of GW(// p )-processes started at Yq = p. 

One important ingredient of the proof is Theorem 3.4 [17] due to Grimvall that is the exact 
analogue of Theorem 1.4 without immigration. For convenience of notation we re-state it as a 
lemma. 

Lemma 4.1 (Theorem 3.4 [17]) The three following assertions are equivalent 

(a) For any t > 0, 

where the process (Z% ; t > 0) is a stochastically continuous process such that 

Vt > , P(Z t > 0) > and P(Z t < oo) = 1. 

(b) There exists a spectrally positive Levy process X = (X t ;t > 0) with exponent ip satisfying (5) 
such that the following convergence 

Pp( Up(X ie -) (45) 

\ p J p~*°° 

holds weakly in R. 

(c) There exists a conservative stochastically continuous CSBP(ip) denoted by Y = (Yt;t > 0) 
started at Yq = 1 such that the convergence 

holds weakly in D([0, oo),R). 

Remark 4.1 Theorem 3.4 [17] is stated with a different scaling: we refer to the proof Theorem 
2.1.1 [11] to derive Lemma 4.1 from Theorem 3.4 [17]. 

Since obviously Theorem 1.4 (Hi) Theorem 1.4(«), we only have to prove (i) (ii) and 
(ii) (Hi). 

Proof of (i) (H): For any t, A G [0, oo) and any p > 1, we set 

u p (t,\) = -p\og(g^ pt] (e-^)) , ^(A) = - 7p log(/W(e- A /P)) 
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b p (t,X) = / ds <p p (u p (s,X)) , d p (t,X) = -logE 

Jo 

and d(t, A) = - logE[exp(-AZ t *)]. First deduce from (1) 

\px 



e K [7p«] 



Up(t,A) + 6p(t,A) = dp(t,A). (47) 



The convergence of Theorem 1.4 (i) combined with Dini's theorem implies the following asser- 
tions: 

• For any A > 0, d(0, A) = xX and d(-, A) is continuous on [0, oo). 

• For any t, A > 0, c/(t, A) G (0, oo) and lim^_> d(t, A) = 0. 

• For any t > 0, cL(i, •) — ► c/(i, •) uniformly on every compact subsets of [0, oo). 

Let T be any denumerable dense subset of (0, oo) and let E be any infinite subset of N. 
By use of Helly's selection theorem combined with Cantor's diagonal procedure, we can find an 
increasing sequence (pk',k > 1) of elements of E, a set of measures (m t ,i G T) on [0, oo) and a 
measure n on [0, oo) such that m t ([0, oo)) < 1, t G T, n([0, oo)) < 1 and such that for all t G T, 

Vr G [0,oo) s.t. m t ({r}) = , lim P^Y?* (1 < r) = m t ([0,r]) , (48) 

Vr G [0,oo) s.t. n({r}) = , lim P^W™ < r) = n([0,r]). (49) 
Define for any A > and any t G T, 

u{t, X) = -log f e- Xy m t (dy) and p(X) = -log [ e~ A %(*/), 

J[0,oo) J[0,oo) 

with the convention — log(0) = oo so that m t = iff u(i, A) = oo for a certain A > 0. By Dini's 
theorem and standard monotonicity arguments, we deduce from (48) and (49) that for any t G T 
the following convergences hold 

Upkitr) , — ► «(*,•) and ^ — ► <p (50) 

as [0, oo]-valued functions uniformly on every compact subsets of the open interval (0, oo). 
We consider two cases: i/O and x = and we first suppose i/O. By (47), we get 

^«p(t,A) < dp(t,A). 

Then we pass to the limit along (£>&;&: > 1) to show that u(t, A) < x _1 d(t,A) < oo , i G T , 
A > 0. Thus m t 7^ for any t G T and it makes sense to define the function b on TU {0} x [0, oo) 

by 

A) = d(t, A) — xu(t, A) if i G T, A > 
and 6(0, A) = 0, A > 0. Deduce from (50) that for any t G T 

&**(*.•)—► KV), (51) 

fe^oo 
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uniformly on every compact subsets of the open interval (0,oo). We first prove the following 
claim. 

Claim 1: There exists to 6 T such that m to ^ So 

Proof of Claim 1: Suppose that u(t, A) = , for any t £ T and any A > 0. Denote by 
q p € [0,1] the smallest solution in [0,1] of g^ p \z) = z. Observe that t — > g^^(z) is non- 
decreasing for < z < q p and non-increasing for q p < z < 1. Thus, for any p > 1 and any A > 0, 
Up(-,X) is monotone. Since T is dense, then a standard monotonicity argument implies that 

u Pk (t,X) — >0 ,t>0, A>0. (52) 

k— >oo 

We take u(-, A) = 0, A > 0, then we also get 

b(t, A) = d(t, A) , i > 0, A > (53) 

and since the 6 Pfc (-,A)'s are non-decreasing and d(-,\) is continuous, (51) holds for any t > 0. 
Now, set s p = 7 p ~ 1 ([7p(i + s)] — [j P t]) and use the Markov property for Y* ,p at time [y p t] to 

get 

d p (s + t, A) = dp(t, Mp(sp, A)) + b p (s p , A). (54) 

Since s Pk — > s, since the i> Pfe (-,A)'s and the u Pfc (-,A)'s are monotone and since their limits are 
continuous, we get 

u Pk( s Pk^) — *■ and b Pk (s Pk ,\) — > d(s,\), t>0, A > 0. 

k^oo k^oo 

Use this to pass to the limit in (54) to get 

d( s + 1, A) = d(s, 0) + b(t, A) = d(t, A), t, s > 0, A > 0. 

It then implies that d(s,X) = xA, s > and the process Z* has to be a constant process which 
contradicts the assumptions of Theorem 1.4 (i). ■ 

Claim 2: n ^ 0. 

Proof of Claim 2: Recall from the proof of Claim 1 that for any p > 1 and any A > 0, n p (-, A) 
is monotone. Since the /W's are non-decreasing, we get for any p > 1 and any t > 0, 

■^i^p (A A u p (i, A)) < 6 p (t, A) < (L(t, A) (55) 

(use (47)) for the right member). Let to £ T satisfying Claim 1. Then, for any A > 0, u{to,X) € 
(0, oo). So it makes sense to pass to the limit in (55) along {p^; k > 1) with t = to. We obtain 

t f (A A u(to, A)) < d(to, A) < oo, (56) 

which implies the claim. ■ 

Claim 3: For any A > 0, b(-, A) extends to a non-decreasing continuous function on [0, oo). 
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Provided that Claim 3 holds, Dini's theorem combined with a monotonicity argument implies 
the following convergence : 

b Pk (t,-) — &(*,•) t e [0, oo) (57) 
holds uniformly on every compact subsets of the open interval (0, oo). 

Proof of Claim 3: Fix p > 1 and A > and recall that u p (-,X) is monotone and that <p p is 
non-decreasing. Thus, we get for any < s < t, 

0<b p (t,X)-b p (s,X) < MziM^ ( u ( t>A ) v«p(a,A)) 

< M_M^^ (dp(t)A)vdp(S)A)) 

by (47). Next, by (51) and Claim 2 

0<b(t,X)-b(s,X)<(t-s)tp(x' 1 {d{t,X)Vd{s,X))) , A > 0, s,teT, (58) 

which completes the proof of the claim. ■ 
Thus, it makes sense to extend the definition of u on [0, oo) x [0, oo) by setting 

u(t, X) := x -1 (d(t, A) - b(t, A)) t > 0, A > (59) 

and u(t, 0) = 0, t > 0. Observe that A) < d(t, A) and n(t, A) < A), for any i > and 

any A > 0. If A goes to 0, then u(t, 0+) = b(t, 0+) = for any t > 0. It implies in particular 
that for any t € T m t is a probability distribution. This result combined with (57) also implies 
that the convergence 

««(*,■),— «(*,■), *G[0,oo) (60) 

holds uniformly on every compact subsets of the closed interval [0, oo). Now, set for any t,y € 

[0,oo), 

Y t (k,y)=p- 1 Y [ ^ kt] with Y** = \p k y]. 

As a consequence of (60) and of the continuity theorem for Laplace exponents (see [14] p. 431), 
there exists a family of probability measures (Pt(y,dz);t,y > 0) on [0, oo) such that the dis- 
tribution of Y t (k,y) converges weakly to Pt(y,dz). In particular, Pt(l,dz) = m t (dz) , t G T. 
Since (t,y) — > J Pt(y, dz) exp(— Az) = exp(— yu(t, A)) is continuous for any A > 0, the mapping 
(t,y) — > Pt(y,B) is measurable for any Borel set i? C [0, oo). Moreover, the Markov property 
and the branching property for the y Pfc 's imply that for any t,s,y,y' > 0: 



/ 



P t (y, dy')P s {y', dz) = P t+S {y, dz) and P t (y, dz) * P t (y', dz) = P t {y + y', dz). (61) 



By Theorem 4 [30] of Silverstein (see also the correspondence between spectrally positive Levy 
processes and CSBPs in Theorems 1 and 2 [23]) there exists a spectrally positive Levy process 
X with exponent ip satisfying (5) such that u(t, X) is the unique nonnegative solution of the 
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differential equation (3). Then the (Pt(y,dz);t,y > 0) are the transition kernels of a non-zero 
conservative CSBPf^). The Levy-Khintchine formula implies that ip is of the form 

V>(A) = a X + (3X 2 [ ir(dr)(e- Xr -l + \rl {r<1} ) (62) 

i(0,oo) v J 

with a Gl,/3>0 and f ir(dr)l A r 2 < oo. 

Concerning the immigration exponent, use (56) to get <^(0+) = d(t,0+) = 0. Thus, n is a 
true probability distribution that has to be infinitely divisible on [0, oo) since it is obtained as 
a weak limit of marginals of rescaled random walks, (p is therefore the Laplace exponent of a 
conservative subordinator denoted by W . It has to be of the form 

(p(X) = kA + / p{dr) (l - e~ Ar ) (63) 

J{0,oo) V ' 

with k > and / p(dr)l A r < oo. Deduce from (47), (50) and (60) that 

d{t, X) = xu(t, X) + f dsip(u(s,X)) , X,t>0. (64) 
Jo 

We now need to show uniqueness for the limiting functions u, ip and ip: let u, tp and Ip be 
obtained by repeating the previous procedure from another denumerable dense subset T C [0, oo) 
and another subsequence (pk', k > 1). Thus, we must have 

d(t, X) = xu(t, A) + / ds if (u(s, X)) = xu(t, X) + / ds!p(u(s, A)) , t, A > 0. 
Jo Jo 

Differentiate twice the latter equation at t = to get 

— xt/j + if = —xtp + tp and (xip' — <p')ip = (xtp' — <p')ip. 
Differentiate the first expression and deduce from the second one the following equation: 

(xij/ - <p')(il> - = 0. (65) 



Suppose that xip' = <p' on a non-empty open interval (a, b). Differentiate twice this expression 
to get 

( p(dr)r 3 e- Xr = -x [ Tr(dr)r 3 e~ Xr , A > 0, 

J(0,oo) J(0,oo) 

by (62) and (63). Thus, /> = 7r = 0, (3 = and xao = k, which imply that Z* is constant: This 
contradicts the assumption of Theorem 1.4 (i). Accordingly, xip' and ip' must differ at aj>oint 
and by continuity xip' ^ <p' on a non-empty open interval. Thus by (65), we get ip = ip and 
<p = <p, which implies the desired uniqueness. ■ 

Thus, we have shown in the i / case that there exist a non-constant CSBP(^): Y = 
(Y t ;t > 0) started at Yq = 1 and a subordinator W with Laplace exponent ip both satisfying (5) 
such that for any denumerable dense subset T C (0, oo) and any infinite subset E C N, we can 



27 



T. Duquesne. 



Immigration Levy trees. 



find an increasing sequence (pk\ k > 1) of elements of E (depending on T) that satisfies for any 
t G T and any A > 0: 



lim E 

k^oo 



E 



and 



lim E 



E 



-XWi 



It easy to prove that these limits imply for any t > that the following convergences 



(d) 



Y t and p^WjP 



1p 



hold in distribution in R. Theorem 1.4 (ii) follows from Lemma 4.1 (a) ==> (6). ■ 

It remains to consider the x = case: Let f , s G [0, oo); for any p > 1, we may and will 
choose s p such that |s p — s| < 1/7 P and [7 p (s p + t)] = [y P s] + [7 p i]. Use the Markov property for 
Y*' p at [7 p s] to get 

dp(s p + i, A) = d p (s, u p (t, A)) + b p (t, A). (66) 
Since b p = d p , we get for any t, s,X, K > 0, any t € T and any p > 1: 



dp(s p + t,X) — d p (t, A) > d p (s, K A u p (t, A)). 



(67) 



Observe next that since b p = d p , the d p (-,A)'s and d(-,A) are non-decreasing. Since d(-,X) is 
continuous, Dini's Theorem implies that d p (s p + t, A) tends to d(s + £, A) when p goes to oo. 
Choose t in T and pass to the limit in (67) along the subsequence (pk, k > 0) to get 

d(a + 1, A) - A) > d(s, K A u(t, A)), t G T, s > 0, A > 0. 

If u(t, A) = oo, then let A go to in the previous inequality to get d(s, K) = 0, for every s, K > 0, 
which contradicts the assumption of Theorem 1.4 (i). Thus, u(t, X) / oo and consequently 
m t ytz 0, for any t G T. Then, use (50) to pass to the limit in (66) and to get 



d( s + 1, A) = d(s, u(t, A)) + d(t, A) , A > 0, s > 0, t £ T. 



(68) 



Fix so > 0. Since ~P(Z SQ > 0) > 0, d(so, •) has to be a continuous increasing mapping. Then, it 
admits a continuous increasing inverse denoted by A: A(d(so, X)) = A, A > 0. We extend the 
definition of u on [0, oo) x [0, oo) by setting 

u(t, A) := A (d(s + t, A) - d(t, A)) . 

Then observe that u(t, 0+) = for any t £ T. Thus, m t ([0, oo)) = 1 , t G T. Now recall that for 
any fixed A > and any p > 1, u p (-, A) is monotone. It implies (60) by a standard monotonicity 
argument. Now, use (56) to deduce </?(A) < oo , A > and </?(0+) = 0, which both imply that n 
is a probability measure on [0, oo). Then, use similar arguments to those used in the case 
to complete the proof of Theorem 1.4 (ii). ■ 

Proof of Theorem 1.4 (ii) (in): The weak convergence of finite dimensional marginals 
is a straightforward consequence of Lemma 4.1 (a) ==> (b) combined with a simple computation 
based on (1). So, it remains to prove tightness. To that end, we adapt an argument of Grimvall 
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[17]: Fix a > and denote by (Yn' p (a);n > 0) a GWl(n p ,u p ) -process started at Y *' p (a) = [pa]. 
Denote by Qa\-) the distribution oi p~ l {Y*' p (a) - [ap]). Theorem 2.2' [16] (see also Lemma 
3.6 [17]) asserts that the sequence of the distributions of the processes p _1 Yj*' p j, p > 1 is a tight 
sequence in D([0, oo), R) if the two following conditions are satisfied: 

(d) For any t > 0, lim M -+oo limsup p _ 00 P (sup < s < t p~ l Y*^ g] > M^j = 0. 

(e) For every compact set C C [0, oo), 

{(Q(f))*>, a€C,p>l} 
is a tight family of probability measures on R. 

Proof of (e): Observe that 

(q?») -**(—) "*y • 

Thus, (e) easily follows from Theorem 1.4 (n). ■ 

Proof of (d): Fix t > 0. Let if be any positive real number. Observe that for any p > 1, any 

A, y > and any s € [0, t], we have 



(p^Y^iy) < K) < exp(K-d p (s,X)) 



(\pv\ 
K u p (s, A) 

< exp [k - ^A A u p (t, \)j 

(use (47) and the monotonicity of the ti p (-,A)'s). Now, since Theorem 1.4 (ii) implies Lemma 
4.1 (c), we get inf{A A u p (t, A), p > 1} > 0. Thus, it proves that for any K > 0, there exists 
M(K) > such that 

P (jr^jG/) > > 1/2 , y > M(K) , s e [0,t] , p > 1. (69) 

Now use the Markov property and (69) to get 

p ( p ~ lY bA x ) >K )^\ P ( <^ t p ~ lY i^] (x) > M(K) ) ' p - 1 

and (d) follows from the one-dimensional marginals convergence. I 



29 



T. Duquesne. 



Immigration Levy trees. 



4.2 Proof of Theorem 1.5. 

We now consider a sequence of (sub) critical offspring distributions (fi p ;p > 1). Recall notations 
of Section 1 and consider a sequence (r p ;p > 1) of GWI(/x p , r p )-trees where fi p and r p satisfy 
(13) and (14). For any p > 1, denote by R p = (H p ; k > 0) and by DP = {D p ; k > 0) resp. the 
height process and the random walk associated with the forest f(r p ) containing the "left bushes" 
of T p . We also denote by H*' p and by D' ,p the processes corresponding to the forest /(r*) that 
contains the "right bushes" of r v . Recall from Section 2.2 the notations (Lk(T p );k > 0) and 
(Lfc(r*); k > 0) , p > 1. We denote by S the function such that 

HP = X(HP,DP,L{t p )) and i? p = T,(H'' P , D 9,p , L(t°)), (70) 

that is specified by (20), (21), (22) and (23). 

We use Corollary 2.5.1 [11] asserting that (13) and (14) imply the joint convergence 

(p~ 1dP M % 1hP M ; * > o) ^ « > o) (71) 

holds in distribution in D([0, oo), M?). We also get an analogous convergence for H* ,p and D* ,p 
since they have the same distribution. 

Next, we use Remark 2.3 and a standard argument to deduce from the right limit of (13) 
that the following convergence 

{p- x L bvt] {T p ) , p-'L^r;) ; t > 0) ^ (U u V t ; t > 0), (72) 

holds in distribution in D([0, oo),R 2 ). Thus, by (71), (72) and Skorohod's representation theorem, 
we may assume that the following convergences 

^(^ lD L.i'^ lfl L.i) f >„ = < x ' ff )' 

and 

lhn (p- l L bpt] { Tp ) , P~ l L bpt] {T' p )) t ^ = (U,V) 

hold a.s. in B([0, oo), M 2 ), where (X,H), (X',H') and (U, V) are independent processes and 
where (X,H) and (X',H') have the same distribution. For convenience of notation, we keep 
denoting in the same way the random processes involved in the latter almost sure convergences, 
so we may also assume that (70) and (8) hold. We first prove 

a.s. in B([0, oo), R). To that end, let us introduce for any w € B([0, oo), R), the right-continuous 
inverse of uj 

S x {uj) = inf{s > : w{s) > x} , x G R 
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(with the convention inf0 = oo). Set V(uj) = {x G R : S x -(lu) < S x (uj)}. It is easy to check 
that uj — > S x {u) is continuous in B([0, oo),R) at any w such that x £ V(u>) (see Proposition 2.11, 
Chapter VI [18]). By (9), the process x — ► S X (U) = U~ x has a.s. continuous sample paths. Then 
it implies a.s. 

lim S x (p _1 L [7p .](rp)) = J/" 1 , z e 



p— *oo 

-1 



Since U is a continuous increasing process, standard arguments imply 

(SAp^L^.^Tp)) ; x > 0) {U x l -x> 0) (74) 

a.s. in B([0,oo),R) (see Theorem 2.15, Chapter VI [18]). Let us set 

a p {n) := inf{£; > : L k (T p ) > 1 - inf D p } , p, n > 1. 

jr'<n J 

Since i — > It = inf s e[o,t] X s is a continuous process, the following convergence 

( inf p- 1 ^; t >o) (lf,t>0) 
a.s. holds uniformly on every compact subsets of [0, oo). Thus, by (74) 

{i^(M);t> o) ^ (uz} t ; t > o) (75) 

p — >oo 1 

a.s. uniformly on every compact subsets of [0, oo). Next, we set 

vr p (n) := inf {A; > : k + a p (k) > n} , p, n > 1. 
Then, by (22) and (23), we get 

HI = n - ir p (n) + H p (n) and a p (vr p (n) - 1) < n - ir p (n) < a p (vr p (n)). (76) 

We easily deduce from (75) that (P7p) _1 7r p ([p7 p -]) a.s. converges to the identity map uniformly 
on every compact subsets of [0, oo). Then, (73) follows from (76) and (75). Use similar argu- 
ments to prove the corresponding convergence for H p . Thus, we have proved that the following 
convergence 

holds in distribution in B([0, oo), E 2 ). The joint convergence of the corresponding contour pro- 
cesses is a consequence of (18) and (19): denote by q p the function associated with t p as defined 
in Section 2.2; set y p (s) = {p~l P )~ l q P {p~f p s); by (18) we get for every T > 0, 



sup 

s<T 



lp C P ,Ar P )-^H P inpyp{s) \<^ + ^ -p \Hl +1 -Hl\^ (78) 



It* , . 1 

-W < 

~ 7p 7p n<Tp 7p ' "" ^ '"' 

in probability by (77). On the other hand, we get from (19) 



sup \y p (s) - £| < J- sup H p + - > (79) 

s<T 2 2p7 p fc<Tp 7p P7p P^°° 
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in probability by (77). Using similar arguments, we prove analogue convergences in probability 
for the right contour processes. Then, 

(% l C 2pip t (r p ), 7p 1 C 2pipt ( Tp ) ■ t > 0) ^ (%, H t ; t > 0). (80) 

The convergence of the sequence p" 1 !^*^ ^ is a consequence of Theorem 1.4: we easily see 
that (13) implies (11) by taking W\ = U\ + V\ and thus, 

fp(fc-i) = Yl r p( fc 'i) ' k ^ L 

To get the desired joint convergence of the contour processes with the GWI process, argue exactly 
as in the proof of corollary 2.5.1 [11] p. 63-64. ■ 
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